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ABSTRACT 

Magnetic fields are usually observed in the quiet Sun as small-scale elements 
that cover the entire solar surface (the ‘salt and pepper’ patterns in line-of-sight 
magnetograms). By using 3D radiative MHD numerical simulations we hnd that 
these helds result from a local dynamo action in the top layers of the convection 
zone, where extremely weak ‘seed’ magnetic helds (e.g., from a 10“® G) can 
locally grow above the mean equipartition held, to a stronger than 2000 G held 
localized in magnetic structures. Our results reveal that the magnetic hux is 
predominantly generated in regions of small-scale helical downhows. We hnd 
that the local dynamo action takes place mostly in a shallow, about 500 km deep, 
subsurface layer, from which the generated held is transported into the deeper 
layers by convective downdrafts. We demonstrate that the observed dominance 
of vertical magnetic helds at the photosphere and horizontal helds above the 
photosphere can be explained by small-scale magnetic loops produced by the 
dynamo. Such small-scale loops play an important role in the structure and 
dynamics of the solar atmosphere and that their detection in observations is 
critical for understanding the local dynamo action on the Sun. 

Subject headings: Sun: photosphere, chromosphere, magnetic helds; Methods: 
numerical; MHD, plasmas, dynamo, turbulence 


Introduction 


The origin of magnetic held generation is a key problem for understanding solar vari¬ 
ability across a wide range of scales. Modern high-resolution observations of global magnetic 
helds by the Helioseismic and Magnetic Imager (HMI) on NASA’s Solar Dynamics Observa¬ 


tory (SDO) (Scherrer et ah 2012), as well as investigations of small-scale magnetic helds in 
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selected areas using the Solar Optical Telescope (SOT) on Hinode (Tsuneta et al. 2008) and 


the balloon-borne SUNRISE observatory (Solanki et ah 2010), have revealed complicated 


dynamics of magnetic fields on the solar surface, and their tendency to self-organize into 
compact magnetic structures. 

Traditionally, the solar dynamo problem is divided according to two time scales: the 
global dynamo, operating on the scale of the 22-year solar cycle and controlling the global 


toroidal and poloidal helds (e.g. Choudhuri et al.||1995 Dikpati & Charbonneau||1999 

Bran- 

denburg et al. 2012 

, and the local dynamo (e.g. Nordlund et al. 

1992 Cattaneo 

1999 

Cattaneo et al.||2003 

Vogler & Schussler||2007 Kitiashvili et al. 2013a 

Rempel||2014), 

which 


operates on the time scales of granulation and super-granulation from a few minutes to a 


few hours and is believed to be responsible for the Sun’s ‘magnetic carpet’ (Schrijver et ah 


1998 Schrijver fc Title]|2002 ). However, the energetic and dynamic connections between the 
two dynamos are not clear. 


In the global dynamo, the kinetic helicity is produced by the action of the Coriolis 
force on convective turbulence (for recent formulations of the global dynamo theory, see 
e.g. Brandenburg et al. 2005 Pipin 2008 Pipin & Kosovichev 2014). However, in the local 


dynamo problem, without global rotation, the mean helicity integrated over a volume is 
zero. This has caused a debate about the existence of a ‘fluctuation dynamo’ driven by 
non-helical forced turbulence. Recent numerical simulations convincingly demonstrate the 
existence of such a dynamo (Cattaneo 1999 Schekochihin et al. 2004b Iskakov et al. 2007). 
The exact physical mechanism of the fluctuating dynamo is not yet understood, and current 
theoretical interpretations include analytical Kazantsev ( |1968 )-type models, models with 
turbulent anisotropy induced by a mean velocity shear (“shear-current” effect, Rogachevskii 
& Kleeorin 2003), and other ideas (for a review, see [Schekochihin et al. 2007). 


It is important to emphasize that the local dynamo mechanism in realistic solar con¬ 
ditions is different from both these cases. While the mean kinetic helicity on local dynamo 
scales is very small (zero in our model), strong stratification and highly anisotropic near¬ 
surface turbulence lead to strong local helical and shearing flows, particularly in downdrafts 
located in the intergranular lanes. Our simulations show that these flows are the primary 
source of the local dynamo in quiet-Sun regions. 


Observational techniques are now able to resolve flows and magnetic fields associated 


with magnetic flux emergence events on very small subgranular scales (e.g. Centeno et al. 

2007 Orozco Suarez et al.|2008 

Martinez Gonzalez & Bellot Rubio[2009 Thornton & Parnell 

2011 Vargas Dominguez et al. 

2014), which probably reflect the local dynamo at work in 


shallow subsurface layers. In addition, resolving magnetic features though observations on 
smaller and smaller scales has renewed interest in the small-scale magnetic field generation 
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problem and raised a question about the existence of unresolved (or ‘hidden’) magnetic flux 


in the quiet Sun (e.g., Stenflo 1982, 2012 Trujillo Bueno et ah 2004; Shchukina & Trujillo 


Bueno 2011). 


Theories of small-scale (or turbulent) dynamo processes have been intensively studied 
in the past. In particular, the similarity between the MHD induction and vorticity equations 


(Batchelor 1950) initiated discussions about magnetic field amplification due to turbulent 


flows. By analogy with vortices when stretching makes them stronger, twisting and stretching 
of magnetic held lines can amplify the held strength. Thus, repeating stretching, twisting, 
and folding of magnetic held lines can provide a local increase in the magnetic hux and its 


recirculation in a dynamo process (Vainshtein & Ruzmaikin 1972). 


Previous studies were mostly done with direct numerical simulations (DNS) performed 


for artihcially forced turbulent hows (e.g., Meneguzzi & Pouquet 1989 Schekochihin et ah 


2004b 

2005 

) and also for convectively driven hows for a wide range of parameters (e.g.. 

Nordlund et al. 

1992 

Gattaneo|1999 

Brandenburg et al. 32012 

). In addition, recent ‘realistic’- 


type radiative MHD simulations have reproduced the solar surface conditions with a high 
degree of realism, and demonstrated that magnetic helds can be quickly amplihed by local 
dynamos in the upper convection zone from a very weak (~ 10“^ G) seed held, leading to 
more than 1 kG magnetic elements (Vogler fc Schussler|2007 ||Kitiashvili et al.j2013a Rempel 


2014). 


In this paper we present new results of 3D radiative MHD simulations of the local 
solar dynamo for various seed held strengths, from 10“® to 10“^ G (in 5 diherent simulation 
runs), and investigate the development and properties of the dynamo process in the realistic 
regime of solar magnetoconvection, when the turbulent magnetic Prandtl number is not 
prescribed a priori as a constant. Our goal is to investigate the primary physical mechanism 
of the magnetic held generation by local dynamo in the realistic solar conditions. A detailed 
statistical study and comparison are planned for a future paper. 


Numerical simulations of local dynamos on the Sun have two specihc features that must 
be taken into account: 1) the solar turbulence is driven by near-surface convective hows in a 
highly stratihed medium; 2) the hydrodynamic and magnetic Reynolds numbers are so high 
that direct numerical simulations (DNS) with fully resolved turbulent scales are not possible. 

The hrst feature requires accurate modeling of solar conditions including radiative trans¬ 
fer with realistic opacities and a realistic equation of state. This means that direct numerical 
simulations of hows with artihcially forced turbulence have limited applicability to the so¬ 
lar dynamo problem. For our simulations, we use a radiative MHD code (“SolarBox”) 
specihcally developed for modeling solar conditions with high degree of realism. The basic 
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formulation is presented in Section 2.1. 


The second feature of solar simulations requires implementation of sub-grid scale (SGS) 
turbulence models. Previous solar magnetoconvection simulations (e.g. Stein et al. 2003 


Vogler fc Schussl^|2007 Rempel||2014 ) used either artihcial or numerical viscosity to model, 
explicitly or implicitly, dissipation due to turbulence at unresolved scales. In this work, 
we use a Large-Eddy Simulation (LES) approach implemented in the SolarBox code. This 
approach is widely used in engineering applications and tested with laboratory experiments 
and is still a topic of extensive study. The SGS models adopted in our solar simulations are 
developed and tested by using DNS MHD simulations including a fluctuating dynamo driven 
by forced turbulence. This approach provides a verihcation of the LES models. Some initial 


results of this project are described by Balarac et al. (2010), and the LES models which 


currently provide the best performance (Balarac, private communication) are described in 
Section 2.1. 


Perhaps the most signihcant implication of the LES approach for dynamo modeling is 
that in our simulations the characteristic turbulent dissipation parameters are not prescribed 
as constants, but modeled according to the local flow and magnetic held properties. This 
means that the critical non-dimensional parameters, such as the turbulent magnetic Prandtl 
number, are not prescribed but change with local conditions. 


In Section 3, for the hrst time, we describe the physical mechanism of individual mag¬ 
netic held generation events in terms of local plasma dynamics. We hnd a strong correlation 
between such events and local helical hows and vortex tubes. It is commonly believed that 
the how helicity induced by rotation is responsible for the generation of the solar global 
(poloidal) magnetic held. No rotation is included in our simulations, but, because of the 
strong stratihcation of the turbulent surface convection, vortex tubes with helical motion 
are generated, mostly in high-speed (near sonic) downdrafts in the intergranular lanes. The 
mechanism of vortex tube generation, dynamics, and associated atmospheric ehects were 
described in our previous papers ( |Kitiashvih et ah 2012a|[b 2013a|[b ). We show here that 
the local turbulent how helicity in the solar magnetoconvection regime plays a key role in 
the small-scale dynamo process. 


In Section 4, we present our results on the statistical relationship between dynamo prop¬ 
erties and the local turbulent magnetic Prandtl number. Such investigations have not been 
performed before. Also, previous investigations were limited to descriptions of stationary 
dynamo states. 

In Section 5, we investigate the ehects of dynamo-generated magnetic helds on the tur¬ 
bulent properties of solar magnetoconvection (back reaction of magnetic helds). A surprising 
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result is that the resulting kinetic energy spectrum is closer to the Kolmogorov law than the 
initial spectrum without magnetic held. We discuss the evolution of the kinetic and magnetic 
energy spectra from seed-held initiation to the dynamo stationary state and also the varia¬ 
tion of the spectra with depth. This allows us to investigate the process of energy transfer 
in the dynamo mechanism. 

In Section 6, we describe the development of local bipolar magnetic structures that 
probably play an important role in the observed ‘salt-and-pepper’ structure of the quiet- 
Sun magnetograms. Another signihcant new result is that the local dynamo can produce 
relatively long-lived compact magnetic structures. Their formation is related to the amplih- 
cation of shearing and twisting hows caused by growing electric currents, and thus magnetic 
forces, around magnetic elements. 


In Section 7, we discuss the topological properties of dynamo-generated magnetic helds 
in the solar atmosphere and show that they form a canopy-like structure consisting of small- 
scale magnetic loops. Such a topology results in a rapid change with height in the solar 
atmosphere of the ratio between the vertical and horizontal held components and can explain 
the discrepancies in the observed ratios (e.g. Danilovic et ah 2010). The mean strength of 


the dynamo-generated magnetic helds corresponds to the observed values from Hinode (Lites 


et al. 2008). In Section 8, we summarize the principal results and discuss their implications 


for our understanding of solar magnetism, as well as future plans. 


2. Computational setup 


2.1. 3D radiative MHD ‘SolarBox’ code 


To perform the local dynamo simulations we use the 3D radiative MHD code ‘SolarBox’, 
developed at the NASA Ames Research Center by A. Wray and N. N. Mansour. The code is 
based on a LES (Large Eddy Simulation) formulation for compressible how and includes a 
fully coupled radiation solver, in which local thermodynamic equilibrium is assumed. Radia¬ 
tive transfer between huid elements is calculated using a 3D multi-spectral-bin method with 
long characteristics. For the initial conditions we use a standard solar model of the interior 
structure and lower atmosphere. The code has been carefully tested and was previously used 
for studying the excitation of solar acoustic oscillations by turbulent convection in the upper 
convection zone ( |Jacoutot et al. 2008bpi ) and several other problems of solar dynamics and 
MHD (e.g. Kitiashvili et al. 2010H 2011, 2013b). 


We solve the grid-cell averaged equations of the conservation of mass ([^, momentum 
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(|^, energy (|^, and magnetic flux (|^: 


dpUi 


dp . . 

^ + {pUi)^i = 0 , 

+ {pUiUj + {Pij + pTij)) ■ = -pet). 


dE 


EUi + (Pij + pTij)Uj - {k + K^)T^i + 


c \2 1 

dvr / cr + cr* 




rad 


dB, 

dt 


+ 


UjBi — UiBj -; ;-— [Bi i — Bii) 

/ ^ 47r(cT + a*) ^ ^ 


= 0 , 


( 1 ) 

( 2 ) 

= 0, (3) 

( 4 ) 




where p is the averaged mass density, Ui is the Favre-averaged (density-weighted) velocity, 
Bi is the magnetic held, and E is the averaged total energy density E = \pUiUi + pe + pef) + 
^BiBi, where 0 is the gravitational potential and e is the Favre-averaged internal energy 
density per unit mass. is the radiative flux, which is calculated by solving the radiative 
transfer equation, and Pij is the averaged stress tensor P^j = {p + \puk,k + ^B^Bk) Sij — 
p{uij + Ujj) — -^BiBj^ where p is the dynamic viscosity. Gas pressure p is a function 

^ " is 


of e and p calculated through a tabulated equation of state (Rogers et ah 1996); Ti 




the Reynolds stress, k is the molecular thermal conductivity; k* is the turbulent thermal 
conductivity, a is the molecular electrical conductivity, and a* = /Anptn is the turbulent 
electrical conductivity in Gaussian units. The magnetic diffusivity, is calculated using the 
subgrid-scale (sgs) model ( Theobald et al.||l994 Balarac et al.||2010 ): = aA^|V x B|/y/p, 

where a is a constant (a = 2), A is the grid spacing, and overbars denote the resolved 
held. The turbulent heat conductivity is given in terms of the turbulent thermal Prandtl 
number, PP, as P = pCpP/PP, where A = /p. We assume that PP is unity. 


Gurrently it is impossible to achieve a realistic Reynolds number in direct numerical 
simulations of solar MHD phenomena; therefore modeling of the dynamical properties of 
solar convection is performed through implementation of subgrid-scale LES (Large Eddy 
Simulation) turbulence models. These models create a realistic turbulent Reynolds number 
at which the simulations can be accurate and so provide a better representation of small- 


scale motions (e.g. 

Zhou et ah 2014). Here we used a Smagorinsky eddy-viscosity model 

(Smagorinsky 

1963 

), in 

which the compressible Reynolds stresses were calculated in the 

form ( 

Moin et ah 

1991 

): Tij = -2CsA^\S\{Sij - Uk,kSij/3) 2CcA^\S\‘^6ij/3, where the 


Smagorinsky coefficients are Cs = Cq = 0.001; Sij is the large-scale stress tensor. 


The simulation results presented here were obtained for a computational domain of 
6.4 X 6.4 X 6.2 Mm, which includes 1-Mm of the lower atmosphere. The grid-size is 12.5 km 
in the horizontal direction. Above the solar surface the vertical grid size is constant (12 km) 
but stretched out in the deeper layers. The lateral boundary conditions are periodic. The 
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top boundary is open for mass, momentum, and energy fluxes, and also for the radiation flux. 
The top boundary is handled through characteristic boundary conditions; hence fluxes of all 
quantities can occur through this boundary in either direction. The bottom boundary is 
handled as follows: 1) the vertical velocity component is set to zero; 2) horizontally uniform 
numerical fluxes are created at the bottom boundary to compensate for any net inward or 
outward flux through the top boundary, for mass, momentum, energy, and B-£eld, except 
that the radiative energy flux through the top is not counted in this process; 3) an energy 
flux equal to the stellar luminosity is introduced at the bottom boundary, also uniformly in 
the horizontal dimensions. As a result, all conserved quantities retain constant integrals over 
the computational volume, except that the total energy is allowed to vary due to radiative 
loss and the steady bottom-boundary flux that mimics energy flow from the core. This seems 
to us to be a physically natural treatment of the energy, and it also serves as a check on the 
radiative transfer calculation. We have verihed that, over a sufficiently long time interval, 
the constant power input at the bottom boundary comes to be compensated by the radiative 
flux through the top boundary, resulting in a total energy that oscillates randomly and with 
very small amplitude around a constant value. In summary: 1) mass, momentum, energy, 
and B-held can advect and diffuse through the top boundary; 2) these losses or gains are 
compensated at the bottom boundary; 3) power to satisfy the stellar luminosity is introduced 
at the bottom, and radiation through the top quickly comes to balance this input within a 
few percent and eventually balances it exactly in time average. 


Compared to early magnetoconvection simulations by Stein et al. (2003) and later work 


by Vogler & Schiissler (2007) and Rempel (2014), our simulation domain is significantly 
deeper: 5.2 Mm vs 2.5 Mm of Stein et al. (2003), 0.9 Mm of Vogler & Schiissler (2007), 
and 2.3 Mm of Rempel (2014) in his high-resolution runs. The deeper domain reduces the 
effects of the bottom boundary conditions. The simulation results presented in Section 3 
show that the primary region of magnetic held generation and energy exchange is located in 
the upper 1 Mm layer, but slowly extends into deeper layers. Therefore, it is important that 
the computational domain has a signihcant depth. 


2.2. Initial conditions 

Local dynamo action is a complicated interaction of magnetic helds and highly turbulent 
hows on small scales. In this paper, the dynamo simulation is started by adding a very 
weak seed held to a hydrodynamic simulation of fully developed solar convection. After 
introducing the ‘seed’ held, no additional magnetic hux is introduced into the domain. To 
investigate the ehects of the initial seed-held properties, we consider 5 cases of magnetic held 































initialization (Table [^. In these cases the initial seed magnetic field has various distributions: 
(A) uniform vertical field, {B) checkerboard-like, alternating polarity vertical field patterns, 
and {C — E) a. randomly distributed net of magnetic field with amplitude ±10“^ — 10“® G 
(white noise for all field components). In case B the checkerboard structure has a period 
of magnetic field variations of 100 km, in order to mix opposite-polarity patches in the 
intergranular lanes. The hydrodynamic conditions at the time of magnetic field initialization 
are exactly the same for cases B — E. 

Table 1: Properties of the seed magnetic field. 


Cases 

Magnetic field 
strength, G 

Initial field 
configuration 

A 

10-2 

vertical 

B 

o 

1 

checkerboard 

C 

o 

1 

to 

white noise 

D 

o 

1 

white noise 

E 

10-6 

white noise 


3. Generation of small-scale fields by turbnlent plasma 


Magnetic field amplification in turbulent solar convection can roughly be divided into 


three basic mechanisms: 1) magnetic field concentration due to converging flows (e.g. Nord- 


lund [1983), 2) convective collapse (Parker 1978 Spruit & Zweibel 1979), and 3) dynamo 


processes driven by helical or shearing motions (e.g. Batchelor 1950 Vainshtein & Ruz- 


maikin 1972; Brandenburg 1995). The turbulent nature of the photospheric layers, where 


the intense radiative cooling drives downward convective motions, makes it impossible to sep¬ 
arate different sources of the locally growing magnetic energy (Spruit 1984), and therefore 
the magnetic field amplification can be discussed only in terms of a dominant mechanism. 


To avoid confusion in separating dynamo and non-dynamo field amplification, no ad¬ 
ditional magnetic flux was introduced into the computational domain after the ‘seed’-field 
initialization. This provides an important test for the dynamo action: if there were no dy¬ 
namo process then local magnetic patches would only be formed due to converging flows in 
the intergranular lanes, and then convective downdrafts would transport these patches into 
deeper layers, from which a part of the flux would be recycled by turbulent motions. In 
contrast, with small-scale dynamo action magnetic field patches are continuously formed, 
resulting in a rapid growth of magnetic energy. 
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We note that during very early times of the ‘dynamo’ runs (shorter than one overturn¬ 
ing time) the evolution of the ‘seed’-held magnetic elements behaves similarly to corks in 
distributed turbulent flows: they tend to collect in the intergranular lanes where the con¬ 
vective flows converge. Figure]^ shows the distribution of the vertical velocity (panel a) and 
magnetic held (panel b) in the photosphere shortly after introducing the ‘seed’ magnetic held 
as a checkerboard structure (case B, see Table [^. Because the weak held mostly follows 
the turbulent how, the regular held distribution becomes deformed. During this stage local 
magnetic held amplihcation is primarily caused through compression by converging hows 
in the intergranular lanes. Shortly after this initial phase we see the start of the dynamo 
process: small-scale helical motions in the intergranular lanes drag magnetic held lines and 
stretch them, and can even reverse the initial local polarity of the held. The dynamo process 
becomes stronger with time as turbulent hows on larger scales become involved. The mean 
magnetic energy density becomes saturated after 7-8 solar hours. Snapshots of the vertical 
velocity and magnetic held in the developed state at z = 0 are shown in Fig. d. The 
maximum magnetic held magnitude reaches 2 kG, and the photospheric mean magnitude 
reaches ~ 20 G (Fig. [^). 

The dynamics of the small-scale dynamo is illustrated in the movie included in the sup¬ 
plementary materials, in which the semi-transparent color scale corresponds to the magnetic 
held strength. The color scale is saturated at 500 G, but the held strength can reach more 
than 2000 G at the photosphere. After the seed-held initialization, the hrst signatures of 
the amplihed held (light green-blue dihuse structures) appear above the photosphere due to 
the short dynamical time-scales there. Shortly thereafter the hrst magnetic patches appear 
at the photosphere. The appearance of opposite-polarity patches reveals a ‘salt-and-pepper’ 
distribution of magnetic helds (Fig. [^). The patches evolve rather chaotically, but never¬ 
theless the maximum magnetic held strength grows continuously until it hnally saturates. 
These magnetic elements are highly twisted and show a very dynamic and non-st at ionary 
behavior, and none of them survive for long. 


Because the local dynamo is driven by turbulent motions that cause magnetic held 
amplihcation due to the stretching and twisting of held lines, it is natural to consider the 
relationship between turbulent plasma motions and the local evolution of the magnetic held. 
For instance, a relationship between the dynamo process and vortex motions produced by the 


Goriolis force in a rotating stratihed medium was demonstrated numerically by Nordlund 


et ah (1992). Our simulations reveal small-scale vortex tubes in the intergranular lanes; 
however, no external rotation has been introduced. These vortex tubes are a result of 
overturning convection and Kelvin-Helmholtz type instabilities (Kitiashvili et al. 2012b). 
They play a signihcant role in the small-scale dynamo mechanism. 
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To quantify the apparent close relation between the appearance of magnetic patches and 
local twisting motions observed in our simulations, we calculate the kinetic helicity density, 
hf = v-(Vx v), and investigate its correlation with the time derivative of the magnetic energy 
density. A similar relationship between the kinetic helicity density and local amplification of 


magnetic flux (analogous to the shear-twist-fold dynamo of Vainshtein & Zel’dovich 1972) 


was previously noticed by Brandenburg et ah (1996). In their simulations the vorticity was 


induced by rotating the simulation box. In our case without global rotation, the mechanism 
of vorticity generation due to shearing instabilities and compression has been discussed in 


our previous paper (Kitiashvili et ah 2012a). 


Vertical prohles of the rms velocity and the kinetic helicity density, 77 = v ■ (V x v), 
displayed in Figure show that the strongest turbulent flows occupy a relatively thin 1-Mm 
deep layer just below the solar surface. Despite the similarity of the vertical profiles for 
velocity and helicity, the kinetic helicity density is more concentrated in the upper 1-Mm 
subsurface layer. Our previous studies of helical motions (vortex tubes) in the near-surface 


layers (Kitiashvili et ah 2012b) showed that the lifetime of the vertically oriented vortex 
tubes can be much longer (up to 1 hour) than the granulation lifetime (~ 10 min). The 
vortex tubes are characterized by strong flows reaching sonic speeds, by a sharp decrease of 
the gas pressure, and by intense radiative cooling. The strong horizontal and vertical vortex 
motions are able to attract small magnetic patches, amplify magnetic held by stretching the 
held lines, and, via downdrafts, transport the held into deeper layers, where the held can be 
further amplihed by compression (Fig. |^. 


To investigate the statistical relationship between twisting motions and magnetic held 
generation, we calculate cross-correlations of the squared magnetic held strength, with 
the kinetic helicity density and with the magnitude of the shear stress, [S'! = (2S'jjS'jj)^/^ 
) (Fig. |3- The result is that the distributed magnetic held patches have 
better correlation with the kinetic helicity density (Fig. [^) than with the magnitude of 
shear stress (Fig. [^). These correlation functions decrease with depth more slowly for the 
kinetic helicity density due to transport of the magnetic held by helical downhows (Fig. |^. 
Thus, we conclude that the small-scale dynamo under solar conditions primary develops due 
to local huctuations of the kinetic helicity density, and that shearing huctuations might also 
contribute to the small-scale dynamo. 


(Moin et al. 1991 


Because magnetic held generation and dissipation are highly inhomogeneous and depth- 
dependent, we analyze areas where the magnetic energy grows and dissipates (Fig. [^. We 
hnd that the time derivative of the magnetic energy density is small above the photosphere, 
while at the photosphere the growth rate of the magnetic energy is faster than dissipation. 
In the subsurface layers the rates of growth and dissipation of magnetic energy are very 
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similar. In layers deeper than 1 Mm, the rates of energy variation slowly decrease and are 
characterized by increasing amplitude of the local energy fluctuations. We plan to investigate 
the atmospheric effects in a separate paper. 


4. Turbulent magnetic Prandtl number and generation of small-scale fields 


Magnetic held amplihcation due to swirling and shearing hows has been investigated 
by many authors (see e.g 


Batchelor 

1950 

Vainshtein & Ruzmaikin 

1972 

Nordlund et ah 


1992 Petrovay & Szakaly 1993; Childress & Gilbert 1995 Brandenburg et ah 1996, 2012). 


Stretching, twisting, and folding of magnetic held lines induce growth of the magnetic hux 
under certain plasma conditions, such as a sufficiently high magnetic Reynolds number. The 
magnetic Prandtl number, which is the ratio of the dihusion coefficients of momentum and 


magnetic held, also plays a critical role in dynamo conditions (Nordlund et ah 1992 Vogler 


& Schiissler 2007). Previously it was shown that magnetic held can be easily amplihed for 
large magnetic Prandtl numbers (> 1), (e.g. see 


Brandenburg et ah 1996 Schekochihin et ah 2004b). However, estimates of the magnetic 
Prandtl number for the Sun, based on molecular dihusivities, vary with depth from ~ 10“^ — 


Pao 

1963 

Meneguzzi & Pouquet 

1989 


10 ^ near the base of the convective zone to ~ 10 ^ at the surface (Rieutord & Rincon 


2010). This suggests that the dynamo cannot operate efficiently since, under these conditions, 
magnetic dihusion is much higher than momentum dihusion. Recently, however, it was 


found that dynamos can work for small magnetic Prandtl numbers (e.g. Ponty et ah 2004 


Schekochihin et ah 2005; Iskakov et ah 2007) and that the transition between the regimes 


where magnetic held is either amplihed or dihused by turbulent plasma motions is described 
by a critical magnetic Prandtl number, Pr^, below which the dynamo processes do not 
occur. For theoretical discussions of this problem in the case of incompressible artihcially 


forced turbulence, see papers by Rogachevskii & Kleeorin (1997); Schekochihin et ah (2004a 


2007); Tobias et ah (2011). 


In our numerical simulations, the molecular scales are unresolved, and LES turbulence 
models are used to estimate the plasma dynamics and dihusivities due to motion on sub-grid 
scales. Therefore, in our simulations the magnetic Prandtl number rehects the properties of 
plasma and magnetic helds only on turbulent scales, and is dehned as the ratio of the local 
turbulent plasma and magnetic dihusion coefficients {Prl^ = Figure illustrates 

typical snapshots of the vertical velocity, magnetic held strength, electric current density, 
the time-derivative of magnetic energy density, the turbulent magnetic Prandtl number, and 
the turbulent magnetic dihusivity in a fraction of the simulation domain at the photosphere 
( 2 ; = 0). The simulations reveal concentration of the magnetic held in the intergranular lanes 
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and highly inhomogeneous time-evolution of the magnetic energy on scales smaller than the 
width of the lanes. In particular, we hnd that the magnetic energy grows mostly at the edge of 
granules and at the periphery of turbulent vortex tubes, where shear motions are strongest, 
whereas the decrease of magnetic energy mostly occurs in areas where strong downflows 
transport magnetic field into the deeper layers. This inhomogeneity is also reflected in the 
distributions of turbulent magnetic Prandtl number Pr\^ and magnetic diffusivity rjl^. Small 
values of are mostly associated with areas where shearing flows are present and the 
magnetic diffusivity has its largest values. 

A statistical distribution of the turbulent magnetic Prandtl number {Pt\^ for a sim¬ 
ulation of 1-hour of solar time is shown in Figures and According to our simulation 
results, Prl^ can vary by several orders of magnitude, and most values are in the range of 
0.1 — 1 (Fig. [^). The local growth and decay of the magnetic energy density correlate with 
turbulent Prandtl values PP^ ~ 10“^ — 1 (Fig. |^). The distribution of the magnetic energy 
density rate of change for different ranges of the magnetic Prandtl number (Fig. shows a 
distinct asymmetry between the decay and amplification of magnetic field, in particular, for 
smaller Pr^. This is clear evidence that the local dynamo processes are particularly efficient 
in areas characterized by low magnetic Prandtl number. 


5. Effect of small-scale dynamo action on the tnrbulent properties of solar 

magnetoconvection 


The quiet Sun regions are characterized by relatively weak mean magnetic flux on the 
solar surface, but this flux does affect the turbulent properties of convection on small scales, 
where magnetic held patches can partially suppress plasma motions. For instance, com¬ 
parison of the turbulent spectra of the photosphere, calculated from simulations with and 
without magnetic held, shows that the slope of the inertial range changes from —11/5 in 


the purely hydrodynamic case to —5/3 in the presence of magnetic helds (Kitiashvili et ah 


2013c). Thus, the dynamo-generated magnetic helds locally act on surrounding turbulent 


hows (the so-called ‘back reaction’). In the current simulations, which account for turbulent 
dynamics on sub-grid scales, we are able to capture the complicated interaction and energy 
exchange between the small-scale helds and hows. 


The dynamo process produces magnetic patches with a wide range of scales of magnetic 
held distribution (Fig. [^) and involves the redistribution of energy among diherent scales. 
Therefore, we consider the development of dynamo as an evolution of the turbulent kinetic 
energy density spectra, which show a redistribution of energy through diherent scales during 
the transition from the ‘seed’ held to a saturated magnetoconvective regime (Fig. |^). In 
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particular, we find that the initial development of magnetic patches on the smallest scales 
increases small-scale electric currents and kinetic energy (red curve on panel a). Then the 
energy shift from small to larger scales changes the slope of the spectra so that they approach 
^-5/3 curve). At the same time, the continuing amplihcation of magnetic helds at 

large wave numbers (small scales) causes suppression of the kinetic energy on these scales. 
Thus, the evolution of the kinetic energy-spectra during the development of the dynamo 
process can be interpreted as partly energy transport from smaller to larger scales and 
partly as a transformation of turbulent kinetic energy into magnetic energy. Note that, 
because of the relatively small size of the computational domain and occasional appearance 
of intermittent ‘bursts’ of the local magnetic held generation (similar to previously observed 


the bursts in DNS simulations by Pratt et al. 2013), variations of the energy spectra can 


be due to transient energy redistributions across turbulent scales (e.g., see the spectrum 
shown by blue curve in Fig. [^). Curiously, in the saturation phase (Fig. [^) the kinetic 
energy density spectrum at the photosphere is of the Kolmogorov type whereas the 


results by Rempel (2014) do not reproduce this slope. Also, the kinetic energy density in 
Rempel’s simulations at k=10 Mm“^ is an order of magnitude higher, and at k=100 Mm“^ 
is almost two orders of magnitude higher compared to our simulations. We do not know the 
reason for this difference, but perhaps the higher kinetic energy density is a result of the 
shallower simulation domain (2.4 Mm) in Rempel’s simulations. The characteristic scale of 
the solar granulation at the photosphere is in agreement with both observations and Rempel’s 
simulations (note, the dehnition of the horizontal scale f = kh/27r in his paper). In deeper 
layers the slope of the spectrum in the inertial range decreases in magnitude: the power law 
is closer to k~^^^ at a depth of 1 Mm, and the slope is even smaller in deeper layers. 

Figure [TO] shows time and depth variations of the power spectra of the magnetic energy 
density from 2 Mm below the solar surface up to the photosphere. We found signihcant 
variations of the magnetic energy spectrum during the hrst half-hour after the ‘seed’ held 
initialization, i.e., during the exponential growth phase. Later in time, the shape of the 
spectrum does not change, and the spectral power shows a steady energy increase at all 
scales. The slope of the magnetic energy spectra for small wavenumbers changes after held 
initialization. It reaches a value of 1/3 after ~ 3 — 4 hours in the subsurface layers and a bit 
earlier in the photosphere. Our simulation results show the existence of a k~^ power law in 
the saturation phase only for depths of about of 2 Mm and below, where the distribution of 


magnetic held is more homogeneous. We do not reproduce the theoretical k law (Golit 


syn 1960 Mohatt 1961[ ). This possibly rehects deviations from the theoretical assumptions 
of magnetic held weakness, and thus the absence of magnetic held back-reaction, in previ¬ 


ous models. Note, that Rempel (2014) didn’t see a clear indication of a power law for the 


magnetic energy in his simulations. Figure 16 shows the averaged over one hour vertical 
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distribution of the Poynting flux and the horizontal magnetic flux density transported by 
upflows and downflows. The magnitude of the Poynting flux is similar to the values found 
in Rempel’s simulations for the open boundary/zero field boundary conditions. However, in 
our simulations the Poynting flux is decreasing at depths below 3 Mm (Fig. 16 r) while in the 
Rempel’s simulations it continues to increase. The reason for this difference is unclear, and 
deserves further detailed investigation. Most of the horizontal flux is transported downflows 
into the deeper convection zone (Fig. [T^). The mean flux recirculating by upflows is several 
times smaller than the flux in downflows. 


6. Formation of small-scale bipolar structures 


Understanding the local small-scale magnetic field generation process is critical for 
studying complex solar MHD problems, such as the interaction of magnetic fields and flows 
through different scales. In addition to the statistical properties described in previous sec¬ 
tions, it is important to consider this process in terms of local physical properties. Here we 
present two characteristic cases of local magnetic held amplihcation. 

The hrst case, demonstrated in Fig. shows the evolution of the vertical velocity, the 
vertical magnetic held (color background), the horizontal velocity held (arrows), the electric 
current density (color map), and enstrophy (contour lines) in a selected 400 km x 400 km 
region of the photosphere with 15 sec cadence. 


The development of a local bipolar magnetic structure is associated with strong (~5 km/s) 
swirling motions, on scales from ~ 12 — 25 km to ~ 300 — 400 km, stretching and twisting 
the magnetic held, which is also compressed by converging hows in the downdrafts. Such 
small-scale swirling motions in the intergranular lanes are typically associated with strong 
downhows, ~ 6 — 8 km/s ( [Kitiashvili et ah' 2010| 2011). The appearance of a bipolar mag¬ 
netic structure (with a prominent negative (blue) polarity in Fig. [T^) is also a result of 
swirling hows driven by a vortex tube oriented along the solar surface (indicated by the 
white arrow in Fig. 11 a), similar to one previously described by Steiner et ah (2010). The 
horizontally oriented vortex tube captures magnetic held lines and drags them into the sub¬ 
surface layers. The dynamics of the positive polarity patch is mostly related to the vertically 
oriented helical motions similar to the process described by Kitiashvili et al. (2010, 2011) who 
showed that the process of spontaneous formation of stable magnetic structures consists of 
two basic steps: (1) formation of small-scale hlamentary magnetic structures associated with 
concentrations of vorticity and whirlpool-type motions, and (2) merging of these structures 
due to vortex attraction, caused by converging downdrafts around magnetic concentration 
below the surface. Such complicated how dynamics is caused by interaction of diherently 
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oriented vortex tubes, which, in fact, have orientations varying in space and only locally can 
be regarded as having ‘vertical’ or ’horizontal’ orientation (Kitiashvili et ah 2012b). The 


generated electric current (Fig. 11c) shows a clear correlation with the swirling dynamics of 
convective flows. The strongest electric current density corresponds to areas on the periph¬ 
ery of helical motions or between the vortices, where stretching of the magnetic held lines is 
the strongest. 

Another example illustrates the development of a bipolar magnetic structure due to he¬ 
lical hows and shearing hows along an intergranular lane. Figure [T^ shows a time sequence 
of the vertical velocity, the vertical magnetic held, the magnitude of the electric current 
density, the time-derivative of the vertical magnetic held with overlaid contour lines of en- 
strophy, and the vertical component of the electric current density (dashed curves correspond 
to negative values). An interesting feature of this example is a diherent scenario for the de¬ 
velopment of magnetic elements of opposite polarity (determined by the sign of the vertical 
held component): the positive-polarity magnetic patch evolves following the ‘classical’ sce¬ 
nario of held amplihcation due to helical hows described in the hrst example, whereas the 
negative-polarity patch starts forming in the intergranular lane mostly due to converging 
hows (magnetic collapse). The evolution of the magnetic elements of this bipolar structure 


is accompanied by locally growing electric current (Fig. 12 c, d). Comparison of the velocity. 


pattern, and time derivative of the vertical magnetic held shows that some swirling motions 
induce magnetic held dissipation, probably due to scattering of the held lines, and often the 
held amplihcation takes place on the periphery of swirls or regions with shearing hows. 


The distribution of vertical electric current density (Fig. 13) shows a strong correlation 
with the horizontal velocity of swirling hows in the subsurface layers. In deeper layers, the 
magnitude of the electric current increases, but its distribution becomes dihuse and does 
not show a clear association with the near-surface dynamics. The kinetic helicity density 
patterns show that the scale of the swirling motions increases with depth, from ~ 50 — 60 km 


at the photospheric layer (Fig. 14) up to 120 km at a depth of 300 km below the photosphere 


(panel d). In the deeper layers, the scale of helical motions continues to increase (to larger 
than 150 km), but the distribution of kinetic helicity density becomes complicated and 


consists of opposite-sign helical hows, which however continue swirling together (Fig. 14 
and disappear at a depth of about 500 km below the solar surface (panel /). 


7. Small-scale dynamo and links to low atmosphere layers 

From previous numerical studies it is known that important dynamical and energetic 
links between subsurface turbulent convective hows and the low atmosphere are established 
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through small-scale vortex tubes. In the presence of magnetic held, the vortex tubes represent 
channels of energy exchange between the convective layers and the chromosphere and can 
result in the heating of chromospheric layers. This may be a source of small-scale spicule-like 


eruptions and Alfven waves (Kitiashvili et ah 2012b, 2013b). 


In addition to these effects, our simulations of the local dynamo show that the anisotropy 
between the vertical and horizontal magnetic held components takes place only in the pho¬ 
tosphere and above. The nature of this anisotropy has been the subject of recent debate 


(e.g 

Drozco Suarez et al.|2007i 

kites et al. 

2008 

Danilovic et al. 

2010 

Ishikawa & Tsuneta 

2010 

Steiner & Rezaei||2012 

Stenho||2013 

). The anisotropy changes with height, but there is 


no dependence on the ‘seed’ held. The vertical distribution of the mean unsigned magnetic 


held components (Fig. 15 1 ) shows a slow increase of their strength with depth that is a 


result of compression. In the convection zone, the distribution of the vertical and horizontal 
magnetic helds is statistically similar. This diherence between the vertical and horizontal 
held strengths is small with a weak dominance of the horizontal helds below the solar sur¬ 
face. Indeed, the ratio between the transverse held, dehned as By, = + B'^, and the 

unsigned vertical held is only slightly above \/2 (Fig. 15?), meaning that the strength of all 
the held components is similar. Near the surface layers the vertical magnetic hux becomes 
dominant. Above the photosphere the mean vertical unsigned hux slowly decreases, whereas 
the horizontal held increases and becomes dominant (Fig. IT^. Note that our results are 


diherent from the previous simulations by Schiissler & Vogler (2008) and Rempel (2014) 


who found that the mean horizontal component of magnetic held is larger than the vertical 
component in the whole range, from the deep photosphere to the chromosphere. Our results 
indicate that the horizontal component is dominant only in the upper photosphere and low 
chromosphere, from about 300 km to 700 km. 

This dominance of the vertical or horizontal helds in the diherent layers rehects the 
topological properties of the dynamo-generated magnetic helds, which are characterized by 
held lines organized in small-scale magnetic loops above the photosphere (Fig. mi' Such a 
topological structure resembles the magnetic canopy suggested from observations ( [GiovaneUi 
1980 Jones & GiovaneUi 1982 Schrijver fc Title]|2002 ). According to our simulation results, 
in the intergranular lanes magnetic held is mostly vertical and becomes horizontal above 
the photosphere. The height of these loops is greater when the magnetic held is stronger 
at the loop footpoints. Such a topology of the magnetic held lines was suggested both by 


observations (e.g. kites et al. 2008) and simulations (e.g. Schussler & Vogler 2008 Steiner 


et al. 2008). The horizontal magnetic helds are stronger at the granule edges and in the 


regions where strong turbulent motions are present. This tendency of the distribution of 
transverse magnetic helds was previously found in observations (kites et al. 2008). For a 


detailed comparison of the simulation results with the observational data it is necessary 
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to perform spectro-polarimetric analyses for the simulated data, taking into account the 
instrumental characteristics. We plan this work for a future paper. 


In addition, it is interesting to note that the closest opposite-polarity patches may not 
even be connected by magnetic field lines above the solar surface as simple magnetic loops, 
but instead interact through electric currents above and below the photosphere forming 
very complicated structures. Figure ^ illustrates the topological structure of the electric 
current density streamlines above and below the photosphere (shown as a horizontal semi¬ 
transparent plane). Each streamline is tracked from a point in the region of positive-polarity 
(orange streamlines) and negative polarity (dark blue) patches. The topological structure 
of the electric currents above the photosphere is often characterized by spirals, arcs, and 
large swirls. Below the solar surface such a topology can represent highly turbulent flows, 
as in the case of the positive patch (orange streamlines, in Fig. 18), or as a very regular 


spiral structure, as in the case of the negative patch (blue lines). For instance, in Figure 18 


a current streamline originating in the positive polarity patch (orange) is strongly twisted 
around the negative patch. 


8. Discussion and conclusion 

In this paper, we addressed the problem of small-scale (local) dynamos that are probably 
responsible for a significant fraction of the background magnetic fields on the solar surface. To 
investigate this problem we used a 3D radiative MHD code and performed several simulation 
runs for different strengths (10“® to 10“^ G) and spatial distributions of the initial seed field 
(Table [^. After injection of the seed field, no magnetic flux was added to or removed from 
the computational domain.We hnd that, in all simulation cases, magnetic held generation is 
qualitatively different only during the hrst few minutes after initialization. This transient 
rehects differences in the initial seed helds. However, after about 4 hours of solar time 
the statistical distributions of the generated helds are similar. Therefore, in this paper we 
mostly focused on analysis of case E, in which the ‘seed’ held was distributed as white noise 
with amplitude ±10“® G. Our simulations show that, due to turbulent dynamo action, the 
magnetic held can be locally magnihed above the equipartition strength (~ 600 G), reaching 
more than 2000 G in the photosphere. 

Our simulation results show that magnetic held amplihcation is driven by converging 
hows into the intergranular lanes, shearing hows, and helical motions. All these mechanisms 
of magnetic held amplihcation are present in our numerical model and linked to each other. 
Thus, the local dynamo process represents a complicated interplay of multiple mechanisms 
that are difficult to separate from each other and that can contribute diherently in individual 





18 


magnetic field amplification events. 


To model the local dynamo process we use the LES approach, which allows us to resolve 
the essential scales of the turbulent solar plasma and model subgrid scales. Unlike in many 
other dynamo simulations, our numerical setup does not restrict the magnetic turbulent 
Prandtl number to a constant value. For instance. Figure [^illustrates a local distribu¬ 
tion of various properties on the solar surface and shows a correlation between the areas of 
magnetic energy growth and small values of the magnetic Prandtl number. The statistical 
distribution of the turbulent magnetic Prandtl number shows that local variations can be of 
several orders of magnitude, but most of the distribution is in the range: ~ 10“^ — 1 

(Fig-H). 

The kinetic and magnetic energy-density spectra during the development of the dynamo 
process show an intense energy exchange between different scales and between the kinetic 
and magnetic energies (Figs [^ and [^. The energy redistribution causes changes in the slope 
of the power law, which in the developed regime reaches the Kolmogorov law 


The development of dynamo action can be seen in the movie (in the supplementary 
material). Shortly after initialization of the ‘seed’ field, magnetic patches appear in the 
photosphere. They are highly twisted and interact with each other. The magnetic energy in 
these patches increases with time. During the evolution of these magnetic elements strong 
helical downflows transport the magnetic held from the subsurface into deeper layers. Some 
of these patches diffuse; however, new magnetic patches are continuously formed. Our 
simulations show that the magnetic patches can have magnetic helds stronger than 2 kG 
(Fig. 1^). 

Our analysis also shows that the local dynamo works most efficiently in the 1-Mm deep 
subsurface layer, where the turbulent hows are strongest (Figs and [^. We conhrm the role 
of vortical motions of plasma in magnetic held amplihcation due to twisting and stretching 
( Nordlund et al.||1992 Brandenburg et ah 1995, 1996). The magnetic patches generated due 
to shearing and vortical motions are transported into deeper layers by downdrafts. In the 
deeper layers the magnetic held can be further amplihed by compression. 


We presented two characteristic examples in detailed studies: a bipolar magnetic struc¬ 
ture generated by the interaction of vertical and horizontal vortex tubes, and magnetic held 
generation by a combined action of a vertical vortex tube and shearing hows. These ex¬ 
amples illustrate three important properties of the process of magnetic held amplihcation. 
The hrst is a strong coupling of various magnetic held amplihcation mechanisms, which 
leads to increased local magnetic energy. The second key property of local dynamo action 
is the multi-scale nature of this process, which evolves in turbulent hows from the smallest 
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resolved scales to the granular scales (Figs 13 and 14). The third property is the complex 
topological and dynamical structure of the process, illustrated in Figure 18, and which is 
reflected in the interaction of individual magnetic patches and surrounding magnetic helds 
and in self-organization into a magnetic network. In this paper we described just some of 
these properties of magnetic held evolution; further investigation is required. The primary 
topology of the dynamo-generated magnetic held is represented by compact magnetic loops 
appearing as bipolar structures in the intergranular lanes (Fig. [^; these loops reach higher 
levels for stronger magnetic held concentrations at the footpoints of the loops. 


Formation of magnetic loops in the solar atmosphere rehects the observed height- 
dependent anisotropy of vertical and transverse small-scale magnetic helds (Fig. 15). We 
found an equipartition of the held component in the layers deeper than 2-Mm below the so¬ 
lar surface, but above the solar surface the transverse magnetic helds are dominant because 
of compact loop topology. This variation of the vertical and horizontal held anisotropy can 
explain discrepancies among diherent observations (e.g Orozco Suarez et ah]|2007 kites et ah 


2008 Danilovic et al. 2010 Ishikawa & Tsuneta 2010 Stenho 2013) and support previous 


numerical analysis (Schiissler & Vogler 2008; Steiner et al. 2008). 


Investigation of the small-scale dynamo action is a key element for understanding dy¬ 
namical and eruptive processes on the Sun and forms one of the building blocks of the 
global dynamics. In this paper we suggest a new view of the problem through analysis of 
realistic-type radiative MHD simulations of the solar subsurface and chromosphere. We plan 
to continue analysis of these data, make detailed comparisons with observations, and clarify 
the physics of the small-scale dynamo action and organization of the Sun’s magnetic network. 
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Fig. 1.— Deformation of the initial checkerboard-distributed magnetic held (panel b) for 
the case B (Bqz = ±10“^ G) due to the surrounding turbulent convection (panel a) at 
the photosphere for f = 30 sec after held initialization. Black-white patterns correspond 
to opposite-polarity magnetic helds, saturated in this image at ±10“^ G. Panels c) and d) 
show distributions of the vertical velocity and magnetic held respectively in the developed 
state, ~ 4 hours after the held initialization, at z = 0 and correspond to the last frame of 
the 3D movie showing volume rendering of the magnetic held strength (see supplementary 
materials). In panel d) the magnetic held image is saturated at ±100 G. 
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Fig. 2.— a) Statistical distribution of the magnetic field strength (|B|) vs. the turbulent 
magnetic Prandtl number (Pr^). Probability distribution function (PDF) of b) Prandtl 
number (the statistical distribution parameters are: mean 0.143, variance 0.067, skewness 
1.83 and kurtosis 2.011); and c) unsigned magnetic held. 
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Fig. 3.— Vertical profiles: a) rms velocity and b) kinetic tielicity density plotted for 1 hour 
with 5 min cadence, 5 h after held initialization. 
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Fig. 4.— Time-depth evolution of the mean helicity density, < H >=< v ■ (V x v) > (panel 
a), and strongest magnetic held strength in a horizontal plane (panel b) for case E (Table 
1) with an initial seed held of 10“® G. 
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Fig. 5.— Time-averaged distribution of cross-correlation coefficients of the kinetic helicity 
density, H (top panels) and the shear stress magnitude, [S'! (bottom panels) with the squared 
magnetic held, at three different depths: 0, —200, and —400 km. 
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Fig. 6.— Profiles of the time-derivative of the mean magnetic energy density divided by the 
total area occupied by positive and negative values, respectively, A = AxAy, as a function of 
depth, plotted separately for the growth (thick curves) and decay (thin curves). Each curve 
is plotted for data taken every 10 minutes, (case E). The dotted curve shows the derivative 
of the total mean magnetic energy density averaged over time. 
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Fig. 7.— Horizontal snapshot of a section of the computational domain in the photosphere 
for: a) vertical velocity; b) magnetic held strength; c) vertical component of the electric cur¬ 
rent density; d) logarithm of the turbulent magnetic Prandtl number; e) turbulent magnetic 
diffusivity; and /) time-derivative of the magnetic energy density. 
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Fig. 8.— a) Statistical distribution of the vs. Lagrangian derivative DEm/Dt. b) and 
c) histograms of DEm/Dt for the various ranges of the magnetic Prandtl number indicated 
in the panels. 
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Fig. 9.— Turbulent spectra of the kinetic energy density. Panel a: time-evolution of the 
kinetic energy spectra at the photosphere. Each curve corresponds to different moment of 
time: 2, 4, 6, and 8 hours after magnetic field initialization. The kinetic energy spectra are 
multiplied by a factor Panel b\ turbulent spectra of the kinetic energy density at the 
saturation state for different depths. Each spectrum is averaged over 30 min. 
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Fig. 10.— Turbulent spectra for the magnetic energy density (case E). Panels a — c show 
the evolution of magnetic energy density for different depths and times indicated in the 
figure. Panel d shows a comparison of the magnetic energy spectra at the saturation state, 
for different depths. Each spectrum is averaged over 30 min. 
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Fig. 11.— Time-sequence with 15 sec cadence in a zoomed 400 km x400 km region, where 
the magnetic held is generated by swirling turbulent hows: a) the vertical velocity in the 
photospheric layer (z = 0); b) the vertical magnetic held evolution showing the development 
of small-scale magnetic elements with opposite polarity (bipolar magnetic structure); black 
arrows represent the horizontal velocity held; and c) the electric current density (background 
image) and the squared magnitude of vorticity (cu^, contour lines). This example corresponds 
to case A, with the initial 10“^ G seed held. The white arrow points to a horizontal vortex 
tube discussed in the text. 
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Fig. 12.— Evolution of a bipolar magnetic structure in the photospheric layer (z = 0) 
illustrated in a sequence of four images with cadence 90 sec for different parameters: a) 
vertical velocity; b) vertical magnetic held; c) magnitude of electric current density; d) and 
e) show the vertical magnetic held growth rate as a background red-blue image; in panel d) 
contours correspond to enstrophy and in panel e) contour lines show the vertical component 
of electric current (dashed curves for negative values). 
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Fig. 13.— The electric current density (color background) and the horizontal velocity held 
(arrows) at different depths, from the photosphere (panel a) to 500 km below the photosphere 


(panel /) for the same moment of time as the first snapshot in Fig. 12 
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Fig. 

field 


14.— The kinetic helicity density, H (color background), and the horizontal velocity 
(arrows) at different depths from the photosphere (panel a) to 500 km below (panel /) 


for the same moment of time as in the hrst snapshot in Fig. 12 
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Fig. 15.— a) Mean vertical profiles of the unsigned magnetic field components (averaged over 
1 hour): vertical component (red curve) and transverse component (blue). Error bars show 
the standard deviation, b) Ratio of the mean vertical and transverse components of magnetic 
field as function of depth below the photosphere. Dashed line in panel b) corresponds to value 

V2. 
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Fig. 16.— a) Mean vertical Poynting flux as a function of depth; b) mean horizontal 
magnetic fluxes in upflows (think curve) and in downflows (thick curve). All properties are 
averaged over 1 hour. Error bars shows the standard deviation. 
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Fig. 17.— Illustration of the topology of the magnetic held lines above the photosphere in 
the local dynamo simulations. The horizontal plane shows the distribution of the vertical 
magnetic held in the photosphere. Red color corresponds to positive polarity, blue color to 
negative polarity of the vertical magnetic held. The range of held strength is from —800 G 
to 300 G. 
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Fig. 18.— Example of the topological structure of the electric current density below and 
above the photosphere. Streamlines correspond to electric currents originating from the 
positive (orange) and negative (blue) polarity patches. The semi-transparent horizontal 
plane shows the vertical magnetic held distribution in the photosphere, where blue color 
indicates negative polarity and red color positive polarity. The snapshot corresponds to the 


rightmost frames of Figure 12 




